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We study the stability and dynamics of an ultra-cold bosonic gas trapped in a toroidal geometry 
and driven by rotation, in the absence of dissipation. We first delineate, via the Bogoliubov mode 
expansion, the regions of stability and the nature of instabilities of the system for both repulsive and 
attractive interaction strengths. To study the response of the system to variations in the rotation 
rate, we introduce a "disorder" potential, breaking the rotational symmetry. We demonstrate the 
breakdown of adiabaticity as the rotation rate is slowly varied and find forced tunneling between 
the system's eigenstates. The non-adiabaticity is signaled by the appearance of a swallow-tail loop 
in the lowest-energy level, a general sign of hysteresis. Then, we show that this system is in one- 
to-one correspondence with a trapped gas in a double-well potential and thus exhibits macroscopic 
quantum self-trapping. Finally, we show that self-trapping is a direct manifestation of the behavior 
of the lowest-energy level. 



I. INTRODUCTION 

Superfluid flow in a toroidal trap is stabilized by a large 
energy barrier between the current-carrying state and a 
state with lower angular momentum jUH]. However, in 
mesoscopic systems, such as atomic Bose-Einstein con- 
densates, the barrier can be sufficiently small that the 
system can tunnel quantum-mechanically to a state of 
lower angular momentum [3J 0]; furthermore, if the short- 
range interparticle interactions are attractive or very 
weakly repulsive, such a barrier does not exist, and the 
system can transition smoothly from the current-carrying 
state, to, e.g., the non-rotating ground state. Recent ex- 
periments in ultra-cold bosonic systems in toroidal traps, 
stimulated by the possibility of shining new light on the 
stability and decay of supercurrents [5] as well as by pos- 
sible applications in other areas, e.g., interferometry [Ql 
and atomtronics [7] , have seen such current decays [8HIU] . 

The stability of superflow depends on the interparti- 
cle interactions, the rotation rate of the trap, disorder 
in the trapping potential, and temperature. We con- 
sider a gas of interacting bosons at zero temperature in a 
toroidal trap rotating at angular velocity £1 and address 
the question of how the single-vortex condensate with a 
metastable or unstable superflow evolves in the absence 
of dissipation, driven either by varying the rotation rate 
of the trap or by varying the interparticle interaction via 
a Feshbach resonance. 

We consider, throughout this paper, a quasi-one- 
dimensional gas of N bosons in a thin annulus of radius 
R and cross-sectional radius ro <C R at zero temperature. 
The basic physics of the stability can be most simply un- 
derstood by considering just two single-particle levels of 
the annulus, the non-rotating state, |0), and the state 
with azimuthal angular momentum h per particle, |1). 
The Hamiltonian of this system has the familiar Nozieres 
form [IT] 

^2 = + JL{Nl + Nl + ANoNJ (1) 



where m is the particle mass, V — l^r^R is the vol- 
ume of the annulus, g is the strength of contact in- 
teractions, and No and N± are the number of particles 
in |0) and |1) respectively with N = Nq + Ni the to- 
tal number of particles. The state with Ni = N is a 
single-vortex state and that with No = N is the ground 
state. Figure [T] which shows the energy per particle as 
a function of Ni for different values of the interparticle 
interaction strength, illustrates the energy barrier that 
appears between the single-vortex state and the non- 
rotating ground state when gN/V > ti 2 /2mR 2 . With 
weakening interaction strength, the barrier decreases, 
and for gN/V < ti 2 /2mR 2 , it disappears, leading to in- 
stability of the single-vortex state. 

We first delineate the regions of stability and the na- 
ture of instabilities of the full system as functions of the 
external rotation frequency of the trap, for both posi- 
tive and negative interaction strengths. In general, the 
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FIG. 1: (Color online) Energy landscape of the two-level 
model as a function of the number of particles in state 
|1). Note the energy barrier between the single- vortex and 
ground states for gN/V > h 2 /2mR 2 (solid line, in black); for 
gN/V = h 2 /2rnR 2 (dashed line, in red), the slope vanishes 
at Ni = N, while for gN/V < h 2 /2mR 2 (dot-dashed line, 
in blue) no barrier exists, indicating instability of the single- 
vortex state. 
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stability of the flow is manifest in the small-amplitude 
Bogoliubov fluctuations about the current-carrying con- 
densate. Starting from a mean-field condensate, we in- 
clude Bogoliubov fluctuations [T^J and find the eigenen- 
ergies of the quasiparticle excitations. With decreasing 
repulsion or trap rotational frequency, an energetic in- 
stability |13[ 114] can appear in the system via excitations 
that decrease the angular momentum of the system by 
one unit; the system can lower its energy by exciting 
these quasiparticles. Moreover, we find a dynamical in- 
stability for sufficiently attractive interactions, where the 
quasiparticle eigenenergy becomes complex [T31 E] and 
the system is driven exponentially rapidly in time away 
from the initial state. For a system to evolve due to an 
energetic or a dynamical instability, the presence of dissi- 
pation is necessary in order to remove energy and angular 
momentum; in this paper, we do not include dissipative 
effects, but will in a future publication. With a knowl- 
edge of the instabilities, we then study simple ground 
states that encompass the underlying physics, consisting 
of the two lowest-lying single-particle states. (Another 
example of how an instability indicates the presence of a 
lower-energy metastable ground state is given in Ref. [15] 
where we studied a rapidly rotating trapped Bose gas in 
the lowest Landau level with a vortex at the center of the 
trap.) 

For the system to feel the presence of the trap, the 
trapping potential must break the rotational symmetry. 
We describe the coupling of the system to the container 
by an asymmetric "disorder" potential stationary in the 
frame rotating at angular velocity ft. Within mean-field 
theory, we determine the stationary states of the con- 
densate formed from the single-particle states |0) and |1) 
and find that for sufficiently large interaction strengths, 
dependent on the disorder potential, the system exhibits 
a non-adiabatic response [171 118] to variations of the ro- 
tation frequency, even if Q is changed arbitrarily slowly. 
This behavior arises from the presence of multiple min- 
ima in the energy landscape (separated by a maximum 
or saddle-point which represents an unstable mode) and 
is characterized by the appearance of a swallow-tail loop 
in the lowest-lying adiabatic energy level and a fold-over 
in the occupation probability of the corresponding state 
as functions of the rotation frequency (see Fig. [5] below). 
The swallow-tail loop implies that the response of the 
system to external rotation exhibits hysteresis [17] . 

Moreover, we show that the quasi-one-dimensional 
Bose-Einstein condensate in a rotating annulus can be 
mapped onto the problem of a condensate trapped in a 
double-well potential with Josephson tunneling between 
the two wells. Therefore, macroscopic quantum phe- 
nomenon of self-trapping in double wells [HI [20] also 
appears in such rotating Bose gases, where the system 
acquires a non-zero time-averaged population difference 
between the two components. The onset of self-trapping, 
which is a steady-state population imbalance, exactly 
corresponds to the behavior of the energy levels discussed 
above. 



We briefly note related theoretical studies in simi- 
lar toroidally trapped systems: Bose condensates with 
dipolar intcrparticlc interactions [21 which induce an ef- 
fective double-well Josephson junction, leading to self- 
trapping [22] ; Bose-Einstein condensates with a modu- 
lated, spatially-dependent scattering length [23] ; and hol- 
low pipe optical waveguides with an azimuthally modu- 
lated refractive index which generates an effective double- 
well potential configuration [24] , 

In Sec. [n] we discuss the stability regime of the con- 
densate by studying the energies of the Bogoliubov ex- 
citations. We then analyze the energy landscape of the 



two-mode system in Sec. Ill and demonstrate a swallow- 
tail loop in the energy of the ground state. We construct 
a mean-field description of this system in the presence 
of the disorder potential in Sec. |IV| The appearance of 
swallow-tail loops and cusps in the energy levels and their 
relation to extrema in the energy landscape are studied 
in Sec. |IV A] Finally, in Sec. |rVB| we discuss the connec- 
tion of this system to a trapped condensate tunneling in a 
double-well potential, and the corresponding connection 
of self-trapping in the double- well system to the behavior 
of the adiabatic energy levels discussed in the previous 
subsection. 



II. STABILITY OF THE GROUND STATE 

For a sufficiently thin annulus, the radial and axial ex- 
citations are frozen out, and the angle around the ring 
becomes the only effective degree of freedom. The nor- 
malized non-interacting single-particle eigenstates of this 
system are <pi(9) = (r\l) = e lW /^/2nR with eigenenergies 
€i = (hi) 2 /2mR 2 where hi is the angular momentum and 
r = (R, 9) is the position vector. The Hamiltonian in the 
laboratory frame is 
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where aj is the annihilation operator for a particle of an- 
gular momentum hj, and g = Aith 2 a/m is the two-body 
contact interaction strength with a the s-wave scatter- 
ing length. The Hamiltonian in the frame rotating at f2 
(denoted by a prime) can be written as [25] 
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where Hq — h/mR 2 is the characteristic scale of rotation 
in the system. The non-interacting single-particle energy 
levels of W , depicted in Fig. [2] are periodic in Q. At this 
stage we do not include the disorder potential. 

In the laboratory frame, the condensate ip c (8,t) obeys 
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FIG. 2: Single-particle energy levels in the rotating frame, 
measured in units of fiflo, as functions of SI. 



the time-dependent Gross-Pitaevskii (GP) equation 
h 2 d 2 
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To determine the stability of the system, we construct 
the normal modes of the condensate by perturbing 
the system around the stationary solution ip c (6,t) = 
e _1 ^ { / ft if> c (0), where /i is the chemical potential. We 
expand the condensate wave function in terms of time- 
dependent modes with angular momentum v measured 
relative to the condensate by writing 



VKM = e-^/ ft [^(0) + <W,t)] 



(5) 



where 



61,(0, t) = e iS ^ ]T [u» tp v (0) e^ t/h - v* v <pl(fi) e^ h ] 

(6) 

with e v the eigenenergies, S(9) the phase of tp c (0), and 
u v and v v complex numbers to be determined. (Follow- 
ing Fetter's notation [TS], we explicitly take the phase 
of the condensate out of the sum, whereas other authors 
include the exponential factor in the definition of excita- 
tion wave functions 26 j.) From now on, for brevity, we 
measure angular momentum in units of K, time in units 
of i^Q -1 , energy in units of MIq, and define the dimcn- 
sionless parameters -q — mRg/2ir 2 h 2 r 2 = 2aR/irr 2 and 

ti = n/n - 

We focus, in particular, on the lowest-energy single- 
vortex state, with a condensate of N c atoms in the state 
|1), for which the GP equation implies that ^ = \ +r]N c . 
The modes are described by the two coupled equa- 
tions IT21 
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from which we find the eigenenergy 
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Note that for \v 2 + r]N c < 0, these energies are com- 
plex, indicating that the condensate is dynamically un- 
stable nig. 

The oscillations of the condensate can also be pictured 
in second-quantization as quasiparticle excitations of the 
condensate. In the usual second-quantized Bogoliubov 
formalism, the coherence factors are given in terms of 
e„ [HIMj by 
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(10) 



and the excitation energy in the rotating frame of a quasi- 
particle carrying v units of angular momentum relative 
to the condensate becomes 



e'M = 1/(1 - fi) + M 



r]N c 



(11) 



Self-consistency dictates that N c + J2u^o \ v "\ ~ N. 

We now analyze the stability of the ground state in 
terms of the normal modes. For weak interactions, 
\r)N c \ -c t, expansion to first order leads to £y(f2) — 
\v 2 + v(l - Q) + r) N c- Thus, at Q = and for repul- 
sive interactions, only e_i ~ — \ + i]N c is negative; the 




energetically unstable 




FIG. 3: (Color online) Stability phase diagram in the rota- 
tion rate - interaction strength plane, for the v = — 1 normal 
mode of a condensate with one unit of angular momentum 
per particle. Energetic instabilities are caused by excitations 
with negative energy, whereas those with complex energies 
lead to dynamical instabilities. The current experiments in 
Refs. [TO] lie in the stable region. 
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v = — 1 mode is energetically unstable and anomalous, 
indicating that the correct ground state has lower angu- 
lar momentum than the original single-vortex state. For 
attractive interactions, e_2 — i]N c is also negative, and 
the v = —2 mode is anomalous as well. 

The general stability phase diagram of the v = — 1 
mode is shown in Fig. [3] in the interaction strength - 
external rotation frequency plane. In the hashed region 
where T]N C < — ~, the quasiparticle energy is complex. 
Note that the regions of dynamical instability and en- 
ergetic instability are in agreement with the arguments 
in the appendix of Ref. [14]. In a dynamically unsta- 
ble mode, where the eigenenergy is complex, one of the 
two components in Eq. ^ grows exponentially in time 
while the other decays exponentially An unstable mode, 
living around a maximum or a saddle-point in the en- 
ergy landscape, hints at the existence of a stable lower- 
energy state, corresponding to a modified condensate. 
However, a small-amplitude analysis does not, in gen- 
eral, reveal the nature of the new stable state (see, e.g., 
Refs. [TH [IB] where such modified condensates are ex- 
plicitly discussed). The solid black line, the solution of 
e' v = 0, shows the critical values of interaction strength 
and rotation frequency needed for stability. For a non- 
interacting system, the v = — 1 mode becomes stable at 
fl = Interestingly, faster rotations shrink the ener- 
getically unstable region and stabilize this mode even for 
weakly attractive interactions. At Q = 1, the energet- 
ically unstable region completely vanishes, and the gas 
becomes stable for r]N c > — \. As mentioned before, due 
to the absence of dissipation in our model, the energy is 
conserved, and the instabilities (although present) fail to 
change the state of the system into one with lower energy 
and lower angular momentum. 

In the experiment in Ref. [TU] , where N ~ 8 x 10 4 atoms 
of 87 Rb with unit circulation were held in a ring trap of 
radius R ~ 9 /im, we find r]N ~ 1.8 x 10 3 . Also, the ex- 
periment in Ref. JS], with N ~ 8 x 10 4 atoms of 23 Na in a 
ring trap of radius R ~ 20 ^im, has r]N ~ 2.9 x 10 3 . The 
initial states in current experiments [9, 10J are within the 
stable regime discussed here, far from encountering any 
energetical or dynamical instabilities [37]. However, by 
suddenly changing the strength of the interparticle inter- 
action from repulsive to sufficiently attractive via a Fes- 
hbach resonance, thereby bringing the system from the 
stable region into the dynamically unstable region, one 
would be able to investigate experimentally the evolution 
of the system in the presence of a dynamical instability. 

The above stability analysis was done for an initial 
condensate in |1). Due to the periodicity of the single- 
particle energy levels with respect to the external rota- 
tion frequency (see Fig. [2| , we can extend the same argu- 
ments easily to condensates in higher angular momentum 
states. For a condensate in the chemical potential is 
fi = \i 2 + rjN c , and the quasiparticle energies become 
eJXfi) = v(j-fl) + \v\ (l^+rjNc) 1 / 2 . Thus, the anoma- 
lous v = — 1 mode (which connects \j — 1) and \j + 1) to 
the condensate) becomes stable at Cl = j — \ for a non- 
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FIG. 4: (Color online) The energy per particle, in units of 
ftfio, of the states \N,0) (labeled 0) and \0,N) (labeled 1) 
and that of the barrier state |6) (dashed line) as functions 
of the rotation frequency, for r/N = 1/4. An energy loop 
(labeled by Lb, £o, and L\) emerges due to the existence of 
a maximum in the energy landscape. 

interacting system; its regions of stability in the presence 
of interactions, for j — 1 < O < j, are identical to those 
shown in Fig. [3] 

III. TWO-MODE APPROXIMATION 

As discussed above, when the energy of the v = — 1 
mode becomes negative, the system, condensed in |1), 
prefers a ground state with smaller angular momentum 
than h per particle. As shown in Fig. [2] over the entire 
range < O < 1, the lowest-lying single-particle states 
are |0) and |1). The states | — 1) and \2) have, in gen- 
eral, much higher energies and are not mixed in with the 
ground state by weak interactions; only for r/N > 1 is 
the mixing significant at 51 ~ and f2 ~ 1. Hence, for 
sufficiently weak interactions (r/N < 1), we can keep only 
|0) and |1) in the description of the system and work in 
a two-mode approximation with the truncated Hamilto- 
nian in the frame rotating at Q, 

U' 2 = {\- Q) Ni + ir/(A^ 2 + Nl + ANoN^ (12) 

where Nj = a^aj. The eigenstates of T-L' 2 are the Fock 
states 

|iV ,AM = 7 j^^) N °(a\) m W*e) (13) 

where |vac) is the vacuum. 

This system has been studied extensively in Ref. [17j : 
we briefly recap the results here. Similar to the non- 
interacting case, in the ground state for Q < g, all the 
particles are condensed into |0), while for Q > \, they 
are condensed into 1 1) . As one finds by extremizing the 
energy, the spectrum also acquires an energy maximum 
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when fi — | < rjN, corresponding to the state 



|6> 



;N + (fi - 



|)/2r?,|iV- 



(fi-D/27,). (14) 



This state can be seen as the maximum in Fig. [T] for 
fi = and rjN > | ; as TyiV decreases below | , the region 
in for which such a maximal state exists shrinks, as 
depicted by the dashed line in Fig. [2] for i]N = j. This 
state acts as a barrier between the two ground states, 
making it energetically expensive for density modulations 
(vortex-induced phase slips) to drive the system from one 
minimum to the other, even as fi is varied. Note the loop 
structure in Fig. [I] indicated by the lines labeled Lb, L , 
and L\i whose effect on the response of the system to 
changes in fi will be discussed in the next section. 



IV. SYMMETRY BREAKING IN MEAN-FIELD 
THEORY 



The angular momentum per particle of the system 
changes according to 

d t (£) = -id t Jdei)* c {e,t)^ c {6,t) = 2vlm[4c ] (18) 

which vanishes, as it should, in the absence of V and also 
when the phase of c\cq equals or tt. 



A. Swallow-tail loops 

We now ask how the system responds dynamically as 
the external rotation rate is varied. As we show, the 
non-linearity inherent in the GP equation leads to forced 
tunneling between the energy levels of the system in the 
presence of V. To see this behavior, we first construct 
the steady-state solutions of the GP equation in the form 



IV'c(O)) where [i is the chemical potential, 



a function of N; then, Eq. il7« gives 



We now turn to the question of how the condensate re- 
sponds dynamically to changes in its rotation rate. In re- 
alistic experiments, the potentials felt by the particles are 
never fully rotationally invariant (see, e.g., Ref. [TO]); the 
breaking of rotational invariance changes the way single- 
particle states of given angular momentum are mixed. 
In order to couple the system to rotations of the trap, 
we include in the Hamiltonian a small "disorder" poten- 
tial, V{6 — fit), which is stationary in the frame rotating 
at fi. Such a potential favors a coherent superposition 
of states (a single condensate) over a fragmented (Fock) 
state [351 130]; we assume the following variational 
form for the time-dependent two-mode condensate wave 
function 

|Vc(t)) = v / iV[c (t)|0) + c 1 (t)|l)] (15) 

where |co(£)| 2 + |ci(i)| 2 = 1. The time-evolution of the 
condensate wave function in the rotating frame is gov- 
erned by the GP equation 



In the two-mode model, we can, with no loss of gener- 
ality, take l>(9) — 2z/cos# with V real and positive, cor- 
responding to a coupling between the states |0) and |1) 
of the form 2A r ^Re[c^Co] , so that the mean-field Hamil- 
tonian becomes 



N 



= (i-fi)| Cl | 2 + |r?iV(l + 2 | Co | 2 | Cl | 2 )+2z/Re[ C * C o]. 

(16) 

With this coupling, the amplitudes obey the two coupled 
differential equations 



id t c = V N[2 - \c \ 2 ]c + Vc u 

id t d = r)N[2- |ci| 2 ]ci + Vc + (| - 0)ci. 



(17) 
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Cl + Vc . 



The occupation probabilities of |0) and |1) as functions 
of fi are 

|C0 ' - {\-n)+2{r 1 N-^- 1 -^ ■ (20) 



The eigenstates, then, have the form 



ID 



( M \ 


|n> = 


\co\\ 


\-\ci\ ' 




\\ci\ 



(21) 



where |I) denotes the ground state branch and |II) the ex- 
cited state branch, since a phase difference of ir between 
Co and c\ minimizes the coupling energy 2A^z/Re[c^co] 
whereas a phase difference of maximizes it. 

The chemical potential (the "adiabatic energy level" in 
the sense of Refs. P3JHB]) ^ s f° un d from the determinant 
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together with Eqs. (20) 



(22) 

fourth-order 



The result is a 
equation for /i, with two to four real solutions depending 
on the values of r/N/2v = A and Cl. The chemical po- 
tential is simply related to the energy per particle in the 
rotating frame by 



E' 



1 



M 2 |ci| 2 



(23) 



The energy levels corresponding to states |I) and |II), 
in general, exhibit an avoided crossing as a function of 
fi due to the presence of disorder. Since fj, and E' are 
simply related by Eq. (23), the physical content of their 



corresponding plots is identical. To illustrate the physics, 
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FIG. 5: (Color online) Adiabatic energy levels (top panels, in black), measured in units of htto, and occupation probabilities 
of |0) (bottom panels, in blue) as functions of Q for rjN = V (left column), r]N = 2V (middle column), and r/N — ?>V (right 
column), with V = 1/5. In all graphs, the lowest energy branch and the corresponding population are indicated by I, and the 
top energy level and its corresponding population by II. The arrows A and B in the right column indicate the discontinuous 
change in the population of |0) and the forced tunneling of particles to |1) as Q is changed past the folding point. 



we plot the behavior in terms of \i since it is graphically 
clearer. The upper panels of Fig.[5]show the real solutions 
for /j, as functions of Ct for selected values of 77 TV and V. As 
shown in the figure, and as we will prove at the end, for 
A < 1, the two energy levels have an avoided crossing at 
the critical value of the rotation rate O c = |; at A = 1, a 
cusp appears in the lower branch at this frequency; and 
for A > 1, the cusp gives birth to a loop in the lower 



branch. The loop discussed earlier in Sec. Ill in the ab- 
sence of the disorder potential (see Fig. [4]) evolves into the 
present loop as the disorder is turned on. Note that at a 
given rotation frequency, Fig. [4] shows either two or three 
states, while Fig. [5] shows two or four states; the extra 
state arises from mixing of the upper maximum-energy 
state with lower-energy states (not shown in Fig.|4|. 

As seen in the lower panels of Fig. [5] the derivative of 
the occupation, |co| , of |0) with respect to Cl diverges 
as the cusp appears in the lowest energy band, and the 
occupation folds over itself (the characteristic S shape 
seen in fold catastrophes [351 [Ml [35] ) as the loop emerges 
for A > 1. The swallow-tail loop indicates hysteresis [T7] 
and a lack of adiabatic evolution with [TH [TS1 [571 HH 
133] in a condensate in an annulus. 

To see the physics of the swallow-tail loop, imagine 
that we prepare the system, with A > 1, on the lower 
branch at f2 = and very slowly increase f2 to avoid any 
tunneling to the other branch. The system will, then, 
follow this branch adiabatically until the point where the 
branch terminates and folds back on itself (at > |). 
Upon further increase of f2, the system is forced to make 
a discontinuous jump either to the lower part of branch 



I (indicated by the arrow A in Fig. [5]) or to the upper 
branch II (indicated by the arrow B). Similarly, the occu- 
pation probability of |0) adiabatically follows the change 
in f2 until the branch starts to fold over itself, at which 
point a sudden change in the population of that state 
becomes inevitable with further increase of the rotation 
frequency, as indicated by the arrows in Fig. [5j In other 
words, a fraction of the particles in |0) are forced to tun- 
nel to |1). Sweeping fl in the other direction forces a 
similar behavior on the system as well. 

The folded-over section of branch I of the occupation 
probability and the respective top part of the swallow- 
tail loop of the lowest-energy level (in the right column 
of Fig. [5]) are inaccessible through a sweep of f2 and cor- 
respond to unstable states. In direct analogy with the 
barrier state discussed earlier, they indicate the presence 
of more than one minimum in the energy landscape, sep- 
arated by a maximum or a saddle-point [171 135) . As we 
show in the next subsection, the appearance of the cusp 
(along with the corresponding divergence of d\c \ /dn 
at Q c ) and the swallow-tail loop (along with the corre- 
sponding fold-over in the level population) are related to 
the macroscopic quantum phenomenon of self-trapping 
or self-locked population imbalance. 



B. Self-trapping 

Interesting quantum phenomena, including Josephson 
effects analogous to those in superconductors and also 
chaotic dynamical behaviors, arise from the dynamical 
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FIG. 6: (Color online) Population difference z (solid line, in red) and phase difference <fr/n (dashed line, in blue) as functions 
of t for AE — 1/2 and (a) A = 1 below the critical value, (b) 7.75, (c) 8.2365 just below the critical value, and (d) 8.2374 the 
critical value. The initial conditions are z(r = 0) = 0.6 and 4>{t = 0) = 0. Note how the oscillatory behavior of z(r) changes 
from purely harmonic to anharmonic as A increases; finally, z(t) and 4>(t) both become time-independent as A reaches the 
critical value. 



behavior of the macroscopic phase difference between the 
two components of the time-dependent condensate (15). 



To see this connection, we recast the time-evolution equa- 
tions (17) in terms of the population difference z = 
l c oT — l c i| (where — l<z<l) and the phase differ- 
ence <p — ao — oil between the two constituent states, 
where ctj = arg[cj]. Then, we find 



d T z = — yl — 



z z sin i 



d T (j) = AE + Az - 



where we rescale the time to t = 2vt and set 



AE={\- Q)/2V. 



(24) 
(25) 

(26) 



These equations are identical to Eqs. (3a) and (3b) of 
Ref. [19] which describes coherent tunneling between 
two Bose-Einstein condensates in a double-well poten- 
tial; therefore, all the results of that paper directly apply 
to the present system. The two levels, |0) and |1), cor- 
respond to the two wells. In the double-well system, an 
applied DC voltage between the two wells induces tun- 
neling and, therefore, an AC particle current 



l = 2Nvd T z 



(27) 



between the two condensates [30J. Analogously, in a 
toroidal trap, an external rotation induces a population 
transfer between the two levels. 

The Hamiltonian ( 16 ) written in terms of 4> and z is 



given, to within a constant term, by 



^2 

NV 



-AEz 



\Kz 2 + ^l-z 2 cos0 (28) 



and is a conserved quantity. Given that the dynamics is 
Hamiltonian, the quantum analog of the Poincare recur- 
rence theorem holds [36], and, therefore, the system is 
inherently periodic in time. 

The time evolution of the system, calculated numeri- 
cally, is shown in Fig. [6] for different A. With increase 
of A for a given initial population imbalance, z(0), the 
oscillations in z(r) change from purely harmonic to an- 
harmonic and a plateau appears in z{t) [the nearly flat 



part of the curve in Fig. [6^c)] . At a critical value of A, de- 
pendent on z(0), the oscillation period becomes infinite, 
and the population imbalance becomes time-independent 
at z(t) = z s [see Fig.[6jd)]. For A larger than this critical 
value, z{t) oscillates in time entirely above z s (not shown 
in Fig. [6]). Similarly, for fixed AE and A, there exists a 
critical value z c of the initial population difference for 
which the oscillations cease and z(t) becomes constant 
after a finite time (the plateau continues indefinitely). 
This evolution to a state with a non-zero time-averaged 
value of z(t) (independent of AE and for A greater than 
or equal to the critical value discussed above) is the ana- 
log of the phenomenon of self-trapping in the double- well 
system [151120]. 

The condition to develop self-trapping is that the two 
time derivatives, d T z and d T 4>, vanish simultaneously. 
From Eq. (24), this requires <j> = or it (although it 
appears that the singular point z = 1 also makes z(t) 
time- independent, the correct solution is actually time- 



dependent [JO], as can be deduced from Eq. (17), and 
is thus unacceptable). From Eq. (25), the steady-state 
value z s is given in terms of AE and A by 



AE + Az, ± z< 



1 



0. 



(29) 



Once the system reaches this plateau, it must stay there 
forever, since the equations of motion are first-order in 
time. The critical initial population difference z c that 
leads to self-trapping can be found from energy conser- 
vation. For an initial phase difference, </>(0), we find 



AEz c + ±A2 
= AE, 



\Kz 



-y/l=~z 



cos </>(0) 



(30) 



As illustrated in Fig. [6] the stationary self-trapped pop- 
ulation imbalance z s is, in general, non-zero and only 
vanishes for AE = 0. 

The steady-state self-trapped solution is, in fact, re- 
lated to the adiabatic energy levels discussed above. The 
stationary value (j> = tt leads to the ground state (branch 

I in Fig. [5]), whereas tf> = gives the excited state (branch 

II in FigT5]). We focus first on the point Q = Cl c at which 
the cusp and the tip of the swallow-tail loop appear and 
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for which AE = 0; choosing the minus sign in Eq. (29) 
(corresponding to the ground state), we find three solu- 
tions 



= 0,±Vl - A" 



(31) 



For A < 1, two solutions are complex and unphysical; 
however, when A = 1, all three solutions become degen- 
erate at z a = 0; and for A > 1, we have three distinct 
solutions. This change in the number of solutions at the 
critical point A = 1 is another indication of the occurence 
of the catastrophe [34] |35] when a swallow-tail loop ap- 
pears. Moreover, Eq. ([29]) also yields 



dz s 



1 



A(l - A 2 ) 



d\cq\- 

an 



(32) 



for AE — 0. For A = 1, this quantity diverges; for 
A < 1, it is negative; and for A > 1, it is positive. To- 
gether, these two quantities exhibit the exact behavior, 
with varying A, seen in the bottom panels of Fig. [5] at 
fl = O c . In general, for non-zero AE, finding z s as a 
function of fl by solving Eq. ([29]) , we indeed see the be- 
havior for | Co | 2 depicted in Fig. ta] Hence, the ceasing 
of the coherent oscillation between the two components 
of the condensate (as the system becomes self-trapped) 
and the appearance of a cusp or a swallow-tail loop in 
the lowest-lying adiabatic energy level are in one-to-one 
correspondence. This proves our previous statement that 
the critical disorder strength for which a cusp or a loop 
first appears is A = 1 or, in other words, V — ijN/2. 



rotation rate of the trap. Describing the coupling of 
the system to the rotation of the trap by a symmetry- 
breaking disorder perturbation to the original Hamilto- 
nian, we investigate the steady-state and the general dy- 
namics of this system in a two-mode mean-field model. 
We find that in the presence of these disorders, swallow- 
tail loops appear in the lowest-lying energy band for suf- 
ficiently strong interaction strengths and, as the rotation 
rate is varied, force a sudden, non-adiabatic change in 
the state of the system and the population of the two 
constituent states. Finally, we investigate the connection 
of this system with a system of two condensates tunnel- 
ing in a double- well potential and find these two systems 
to have identical dynamics; therefore, the analog of the 
phenomenon of self-trapping also appears in the system 
studied in this paper. We calculate the onset and the 
properties of self-trapping, and show how the steady- 
state self-trapped states are described in terms of the 
energy eigenstates. 

The next step, which we discuss in a future publica- 
tion, is to include the effects of non-zero temperature on 
the detailed time-dependence of a Bose gas rotating in an 
annulus. At finite temperature, the angular momentum 
and energy of the condensate need not be conserved. As 
a result, thermal, as well as quantum, fluctuations can in- 
duce transitions between metastable states, as has been 
observed experimentally [10J . Furthermore, the coupling 
of an unstable condensate to the environment would al- 
low it to evolve to a stable configuration. 
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